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We theoretically study specific schemes for performing a 
fundamental two-qubit quantum gate via controlled atomic 
collisions by switching microscopic potentials. In particular 
we calculate the fidelity of a gate operation for a configura- 
tion where a potential barrier between two atoms is instan- 
taneously removed and restored after a certain time. Pos- 
sible implementations could be based on microtraps created 
by magnetic and electric fields, or potentials induced by laser 
light. 



I. INTRODUCTION 

The creation and manipulation of many-particle en- 
tangled states offers new perspectives for the investiga- 
tion of fundamental questions of quantum mechanics, and 
is the basis of applications such as quantum informa- 
tion processing. Several proposals to implement quantum 
logic have been made including ion-traps 0, cavity 
QED and photons ||, and molecules in the context of 
NMR Q. Very recently, we have identified a new way 
of entangling neutral atoms by using cold controlled col- 
lisions II (see also ||). Neutral atoms are good candi- 
dates for quantum information processing, since they suf- 
fer a comparatively weak dissipative coupling to the en- 
vironment. Techniques to cool and trap atoms by means 
of magnetic and optical potentials have been developed 
in the context of laser cooling and trapping, and Bose- 
Einstein condensation (BEC) 0. In particular the on- 
going development of magnetic microtraps ||] offers an 
interesting new perspective for storing and manipulating 
arrays of atoms || |nj and possible applications in quan- 



tum information [ ll| 

Motivated by these new experimental possibilities we 
will study in this paper specific configurations of atoms 
stored in time dependent microtraps. We will assume 
that two internal states of the atoms \a) and \b) represent 
the logical states |0) and |1) of a qubit, respectively. The 
aim is to implement a fundamental two-qubit quantum 
gate between two atoms with the truth table 



-|0>|0>, 
|0)|1)-|0>|1), 
|1>|0>-|1>|0>, 

11)11)^-11)11), 



by switching the trapping parameters. Eq. (|l|) represents 
a so called phase gate. To realize this transformation 
we will consider state selective switching of the trapping 
potential such that the atoms pick up a phase due to 
collisional interaction only if they are in state |6). 
This can be achieved by raising and lowering a potential 
barrier between the two atoms as shown in Fig. [j]. Ac- 
cording to Fig. [l]a the potential is initially composed of 
two separated wells. Ideally the atoms have been cooled 
to the vibrational ground states of the two wells. At 
time t = the shape of the trapping potential is changed 
for particles in state \b) (dashed line in Fig. while 
the potentials for the atoms in the state \a) remains un- 
changed (solid line in Fig. |l|b). By removing the barrier 
the particles in state \b) start to oscillate and will collide. 



a) 




-x 




b) 




(i) 



FIG. 1. Configuration at times t < 0, t > r (a) and during 
the gate operation (b). The solid (dashed) curves show the 
potentials for particles in the internal state \a) (\b)). 



The "cold" collision represents a coherent interaction de- 
scribed by a pseudo-potential with a strength propor- 
tional to the s-wave scattering length ||. This results 
in a phase shift of the wave function for both atoms in 
the internal state \b). The size of the phase shift can 
be controlled by the number of oscillations a nd th e ef- 
fective collisional interaction strength (see Sec. II A| ). As 
a last step the atoms have to be restored to the mo- 
tional ground state of the trapping potential of Fig. 0a. 
This whole process of switching the potentials can be 
performed either as (i) switching the shape of the po- 
tential instantaneously at times t — and t = r where 



1 



r is a multiple of the oscillation period in the well of 
Fig. 0b (dashed line), or, (ii) deforming the shape of the 
potential between Fig. [l]a and b adiabatically. The aim of 
the present paper is to investigate the gate dynamics for 
the scenario (i), when the switching is instantaneous: In 
particular we are interested in the required physical pa- 
rameters and the corresponding fidelities characterising 
the quality of the phase gate. We will also study the de- 
pendence of the fidelity on the temperature of the atoms. 
The paper is organized as follows. Section [n] describes 
the model and deri ves an expression for the collisional 
phase shift. In Sec. Ill we study the gate dynamics for 
the case of instantaneous switching while in Sec. IV we 
present numerical results for the fidelity. 



II. MODEL 

In the present section we will write down the Hamilto- 
nian for two interacting particles trapped in conservative 
time dependent potentials and derive an expression for 
the collisional phase shift. 



A. Hamiltonian 

The dynamics of atoms in a time-varying, state- 
dependent trapping potential V a (x, t) (where t is time 
and x = (x, y, z) is the three-dimensional coordinate) 
can be described by the Hamiltonian operator 



H 



J2 Jd 3 x¥ a ( X ) 



ae{a,b} 



2m 



V 2 + K(x,t) 



x E/ a/3 (x,x / )$ /J (x')$a(x), (2) 



a,/3€{a,fc} 



where m is the mass of the atoms, ^ Q (x) is a field oper- 
ator for atoms in internal state |a), and U a ff{pc, x') is the 
potential for the interaction between two atoms in states 
| a) and \/3), where a,f3 € {a, b}. We take a trapping 
potential of the form 



V a (x, t) = v a (x, t) + v± (y) + v± (z) , 



(3) 



i.e. we assume the same shape along y and z, which is 
independent of time and of the internal state. 

For cold atoms the dominant collisional interaction is 
the s-wave scattering term described by a contact poten- 
tial of the form 



t/ Q/3 (x,x') = s - (5 3 (x-x'), 



(4) 



where a"" is the s-wave scattering length for the corre- 
sponding internal states. Note that for identical atoms 



in the same internal state s-wave scattering is only pos- 
sible for bosonic atoms (cf. the b-b collision in Fig. [j]b). 
We therefore require, in the following, the field operators 
^/ Q (x) to describe bosonic atoms and to obey the usual 
bosonic commutation relations. 

Furthermore, we assume much stronger confinement 
along the y and z directions than in a;, so that the prob- 
ability of transverse excitations can be neglected. If each 
atom is initially in the ground state \if>±) of the trans- 
verse potentials, it will then remain in that state to a 
good approximation and the corresponding degrees of 
freedom can be integrated out. In this case the dynamics 
becomes effectively one-dimensional and is described by 
the Hamiltonian operator 



h 2 d 2 , ■ 

2^5 + ««■(*>*) 



h x = Y / 

+ Y \ J dx dx' Tpi(x)Tpl(x' 



tp a (x) 



a,f3e{a,b} 



XU a j3(x-x')lpp(x')4> a (x). (5) 



Here ip a {x) is the one-dimensional analogue of ^ Q (x), 
and 

Uapix — x') — I dy dy' dz dz' i7 Q| g(x, x') 



x W±{y)il)±{y')il)x{z)il>±(z')\* 



S(x - x') 



dy\^±(y)\ A 



(6) 



is an effective interaction potential taking into account 
the transverse confinement of the atoms. The ip± are 
the ground-state wavefunctions in the transverse direc- 
tions (having energy hui±/2 each). Their time evolution 
will just contribute an overall phase factor (with a phase 
proportional to u±), irrelevant for the quantities we are 
going to compute. We sec that the effective interaction 
strength can be adjusted by changing the trapping pa- 
rameters. 

Eq. (||) holds for an arbitrary number of atoms. We 
now consider the case of two bosonic atoms 1 and 2, 
with internal states \a)i.2 and |6)i.2- Their evolution is 
governed by the first-quantized Hamiltonian 



n 



Y n <*f) ® \a)i(a\ ® |/0) 3 </3|, 

a,/36{a,b} 



where 



K 



a/3 



-i/O 

n af3 



~Ha0 + u af3, 

H a (pi,Xl,t) +H/3(p 2 ,X2,t), 



Pi 

2m 



(7) 



(8a) 
(8b) 

(8c) 



Here X\ and pi are the position and momentum operator 
for particle i — 1,2 respectively. 
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B. Phase shift due to interaction 

We call |^g(t)) and \i) a p(t)) the two-particle states at 
time t evolved from the same initial state \il> a p(0)) in the 
absence and in the presence of interaction, respectively: 



ihd t \4> a f3(t)) = H a p\lp a p(t)) ■ 

We also define the overlaps 

O(^ a0 ,t) EE (V»a/S(*)|^(0)). 



(9a) 
(9b) 



(10a) 
(10b) 



The condition that both atoms end up at time t = r with 
the same spatial distribution they had at the beginning 
will not be exactly fulfilled in realistic situations. How- 
ever, in order for our scheme to work it is required that 
this is true at least approximately: 



\O(iP a0 ,r)\ 



1 



V a, 0, 



(11) 



i.e. the two-atom final state should differ from the initial 
one just by a phase factor $ Q| a(T) = axg[0($) a p, r)]: 



|VW»( T )) 



-i*a^(T) 



I^Mo)>. 



(12) 



We also assume that the interaction between atoms does 
not induce any significant alteration in the shape of the 
wave functions, i.e. 



|Oo(^a/s,t)|«l Va,j8,t. 

Hence 

|V«/j(*)>«e- < ^W|^ (t))j 
having defined the collisional phase 

(j> a p(t) = axg,[Oo(i/> a p,t)], 



(13) 



(14) 



(15) 



accounting for the contribution of the interaction to the 
total phase $0/3 (r). The rest of the phase comes from the 
motion of the particles in the time-dependent trapping 
potential. From Eqs. ( |TT| ) and (jl^) it follows that 



V a, /3, 



(16) 



which implies, by analogy with Eq. (12) 



l^(r))«e-^-W + ^W]|^ a/s ( )). (17) 

Here the kinematic phase 4> a {j) [<Pp(t)\ is defined as the 
phase that one atom would acquire after evolving for a 
time r in the potential v a [vp] in the absence of the other 
particle. By substituting Eq. (|rij ) into Eq. (14) evaluated 
at t = t, and comparing it with Eq. (|l2|), the collisional 
phase can be reexpressed as 



Mr) « $ a/J (r) - y> a (r) + M^)]- (18) 
By combining Eqs. (|§), (jl|) and (0), we find 

nd t <f> a p(t) « (vd>)K^3 (*)) = Af?^(t), (19) 

which is precisely the result one would expect from per- 
turbation theory. In order for Eq. (13) to hold, the time- 
dependent energy shift defined in Eq. ( |l9| ) has to satisfy 
the condition AE a p(t) <C Tiuj, with hu> the first excita- 
tion energy of the system. Integration of Eq. (19) gives 
a perturbative expression for the collisional phase: 



4>ap{t) 



dt'AE a p(t') 



(20) 



III. GATE OPERATION 

To proceed further, we have to specify the functional 
form of the potential v a (x, t) in Eq. (|^). The two atoms 
are initially trapped along x in two separate harmonic 
wells of frequency luq, centered at zLxq. In order to sim- 
plify the analytic calculations, the confinement in the 
transverse directions is also assumed to be harmonic. At 
t = the barrier between the wells is suddenly removed 
in a selective way for atoms in internal state \b): an atom 
in state \a) feels no change, whereas one in state |6) finds 
itself in a new harmonic potential, centered on x = with 
frequency lu < loq. The atoms are allowed to oscillate for 
some time, and then at t — r the barrier is suddenly 
raised again to trap them at the original positions. Dur- 
ing this process the atoms acquire a kinematic phase due 
to their oscillations within the wells, and also - if they 
collide - an interaction phase due to the collision. Here 
we calculate these phases and consider the appropriate 
switching time r for a quantum gate. In Sect. IV we 
make a quantitative estimate of the gate fidelity. 



A. Switching potential 



We take the potential in Eq. (||) to be explicitly 



v a (x,t) 



Vb(x,t) 



v±(y) 



mujQ 



[9(x)(x - x Q ) 2 + 8{-x){x + x ) 2 ] ; (21a) 



V a (x,t) t<0, t> 



moj 



-x 2 < t < r; 



y , 



(21b) 



(21c) 



as shown in Fig . [l|. As l ong as the single- well ground-state 
width ao = ^Jh/muoQ satisfies a§ <C x\ and there are 
no significant excitations to higher levels of v a (x,t), the 
actual behavior of that potential arou nd th e origin does 
not really matter and we can use Eq. (21a) regardless of 
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the experimental shape of the barrier around x = 0. The 
ground state wavefunctions ip±(x) of the right and left 
well of the potential v a {x,t) are given by 



(22) 



while the ground state wavefunction in the transverse 
directions is given by 



ip±(y) = hr-r e 2R 



V 2tt?i / 



(23) 



By assumption the overlap between the two wavefunc- 
tions i/j+(x) and ip-(x) is negligible since the two parti- 
cles are kept separated from each other in the potential 
v a (x,t). At t = 0, the central barrier between the two 
wells is selectively switched off for state \b). A particle in 
this state will start moving towards the other atom along 
x and an interaction will take place. We shall separately 
study the evolution of the system at t > for each com- 
bination of internal states (a, (3). For operation of the 
quantum gate analyzed here, it is important that Vb(x,t) 
be accurately harmonic while < t < r. 



B. Particles in the same internal state 



1. Initial state 

If both particles are in the same internal state \a), this 
factorizes from the motional degrees of freedom and the 
initial state is 

\ip aa {0)} = -j= ®|a)|a). (24) 

The calculation can be simplified by introducing the cen- 
ter of mass (CM) and relative coordinates for the x— 
motion, thus rewriting 

1paa(Xl,X 2 ,0) = —f= [^-(xi)lp + (x 2 ) + ll>+(xi)lp-(x2)] 



1/>CM(R, 0)^rel(r, 0), 



where 



77 <, = -!, + 1 



(25) 

(26a) 
(26b) 



with M = 2m, /i = m/2, R — [x\ + x 2 ) /2, r = x 2 — x\ 



2. Time evolution 

For t < 0, the particles are stored in the displaced 
wells and no interaction takes place. If both particles 
are in state |a), the potential remains unchanged also for 
t > 0; there is no collision and thus the collisional phase 
4>aa — 0. The state simply picks up the phase due to the 
free evolution: 



IVwW) =e~^'|Vw(0)} 



(27) 



We shall now consider the situation in which both par- 
ticles are in state \b). In this case, after the barrier is 
switched off, the particles start oscillating in the har- 
monic trapping potential. In the absence of interaction 
they would come back to the initial state after an os- 
cillation period T osc = 2ir/tu, having acquired a phase 
4t:luj_/uj because of the transverse confining potential. 
The interaction causes an additional phase to be accu- 
mulated by the wavefunction as the number of oscilla- 
tions increases, and a slight decrease in the oscillation 
frequency, because the atoms acquire a small delay in 
their motion inside the trap as they come out from a col- 
lision. If the latter feature is not too strong, by choosing 
a switching time r i=s 2Nir/u> it should be possible to get 
back the original state plus an interaction phase, that is 
adjusted to ±7r by a proper choice of the trap parameters 
and of the number of collisions occurring during the ac- 
tual gate operation, i.e. for < t < t. We shall therefore 
focus on the dynamics in this time interval. 

In the center of mass-relative coordinate system we get 



K 



hi, 



2M 



Mlo 2 



R 2 + %- + ^r 2 +u bb (r), (28) 
2jU 



flLU 

~2 



where P — pi + p 2 , p = {p\ — p 2 ) /2. If the interaction 
is neglected we can solve the two-particle Schrodinger 
equation for Ham iltonian Eq. (|28|) analytically a s shown 



in Appendix A 1 



It can be seen from Eqs. ( |All - [A7| ) 
that the unperturbed two-atom motion has a period of 
T osc /2 instead of T osc . This happens because the initial 
state, symmetric with respect to the origin, has nonzero 
projection only on the even eigenstates, having energies 
(2n + l/2)fuj: therefore, after a time 7r/w, each compo- 
nent of the wavefunction gets the same constant phase 
exp[i(2n+ l/2)n] = exp(i7r/2). This has a simple physi- 
cal interpretation: if the atoms do not interact, after half 
an oscillation period each particle is at its turning point, 
coinciding with the other atom's starting location; so at 
that time the two atoms have interchanged their posi- 
tions, but since they are indistinguishable this has to be 
regarded as exactly the same motional state they had at 
the beginning (apart from a phase factor). 

When we take into account the interaction between 
particles, the center of mass motion is unaffected but the 
relative motion can no longer be treated analytically. The 
numerical method we use to carry out this calculation is 
outlined in Appendix A 2 a. It is, however, possible to 
take the interaction into account perturbatively as shown 
in the following section. 



4 



3. Perturbative calculation of the phase shift 

Eqs. (Q) and (|2^) combine to yield 

u a p{xi - x 2 ) = 2a'^ l3 huj_ L 5(x 1 - x 2 ). (29) 

When both particles are in state \b), the time-dependent 
energy shift defined in Eq. © can be calculated analyt- 
ically: 

AE bb {t)= [ dRdr\^ C M(R,t)^(r,t)\ 2 u bb (r) (30) 



J*hu± 



nfi 



where Q(t) is defined in Eq. (A.2). The corresponding 
interaction-induced phase shift accumulated after an os- 
cillation period is 



<fibb(T sc) 



4afu± 



xjx 2 uj 2 - a o^o/ 4 ' 



(31) 



which has been evaluated by means of the well-known 
saddle-point approximation. 



C. Particles in different internal states 

1. Initial state 

When the internal states of the atoms are different, 
they no longer factorize as in Eq. ( pi[ ) and the initial 
state is given by 

|Va fc (0)) = [|V-)i|^+>2 ® \a)i\b) 2 + (1 <- 2)] , (32) 

where without loss of generality we assumed that the 
particle in the left (right) well is in internal state \a) 

(!&))■ 



2. Time evolution 

The relevant quantities can again be expressed in terms 
of the projection of the evolved state on the initial one. 
By virtue of symmetry under particle interchange, this 
turns out to be 



0(^„i,,*) = (V-|(^+|e- 



r'Habt 



I^->IVm 



(33) 



Therefore we can restrict our analysis, as in the previ- 
ous case, to the one-dimensional motion, starting from 
the non-symmetrized wavefunction ip-(xi)ip+(x2)- The 
Hamiltonian for < t < r reads 



pi p 2 2 rrujfi 2 ma; 2 2 

+ u ab (xi - x 2 ) 
P 2 p 2 m, o m i i f 

= 2M + | + 2<- -^ &+ 2*S( 1 -^) 



(34) 



where w = y (^ 2 + Wq) /2. Only the left well of v a (xi,t) 
has been considered since the wavefunction remains neg- 
ligible in the region x\ > for t > 0, as it is at t = 0. 
It can be seen from Eq. (Q) that the center of mass no 
longer decouples from the relative motion, unlike in the 
previous symmetrical case. A numerical calculation is 
needed to evaluate the phase shift <j) ab . This is done in 
Appendix A 2b. 



D. Particles at finite temperature 

Up to now we have assumed the particles to be in a 
well known motional state. In realistic experimental sit- 
uations this may not be the case. The temperature T of 
the particles in the trap will be different from and thus 
the initial state of the system with particles in internal 
states a, (3 is given by the density operator 



Pap(T,t = ) cx 



-n a /3(o-)/k B T 



(35) 



This takes the average over different initial excited states, 
with a thermal probability distribution corresponding to 
T. As shown in Appendix |b] the collisional phase accu- 
mulated is independent of the shape of the wave function 
if the particles move at a constant velocity with respect 
to each other and the shape of the one particle wave- 
function does not change during the interaction. This is 
a good approximation for the interaction between parti- 
cles in the same internal state \b). The particles interact 
in the vicinity of the center of the well where their ve- 
locity v w xquj is almost constant and the shape of the 
one particle wavefunction does not change substantially 
as long as the conditions 



a x < x , 



and 



a x 



(36) 



hold, where a is the width of the one particle wavefunc- 
tion when the particles cross the center of the trap and 
a x = y/h/muj. Therefore the collisional phase </>&fc(T osc ) 
is almost independent of the temperature T as long as 
mainly excitations fulfilling conditions Eq. ( |36"|) are popu- 
lated. Note that we are neglecting transverse excitations. 
If all three motional degrees of freedom are character- 
ized by the same temperature T, this is realistic as long 
as the condition fc^T < fru± is satisfied. However, in 
principle it is also possible to cool the transverse motion 
separately, allowing a higher temperature along x. Of 
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course this would require that the rethermalization time 
is much larger than the experimental time scale. This 
lack of sensitivity to temperature applies quite generally, 
for example to atoms interacting in an optical lattice as 
discussed in H , provided that the velocity at which the 
atoms are made to interact (in that case the velocity of 
lattice movements) is kept constant during the interac- 
tion. 



IV. A PHYSICAL IMPLEMENTATION 

We now consider the implementation of a switching po- 
tential by means of static electric and magnetic trapping 
forces. We first discuss the possibility of obtaining the 
desired state dependence by means of devices which are 
experimentally available ||[ll], when the present mag- 
netic devices can be combined with nanofabricated elec- 
trodes. Then we compute the performance of a quantum 
gate for realistic trapping parameters. 



A. Microscopic electromagnetic trapping potential 

The interaction between the magnetic dipole moment 
of an atom in some hyperfine state \F, mp) and an ex- 
ternal static magnetic field B entails an energy U magn ~ 
<7f/^bw_f|B| depending on the atomic internal state via 
the quantum number tuf (here /is is the Bohr magne- 
ton and <7f is the Lande factor). The Stark shift in- 
duced on an atom by an electric field E gives an energy 
(independent on the hyperfine sublevel) U e \ ~ ia i|E| 2 , 
where a c \ is the atomic polarizability. The interplay be- 
tween these two effects can be exploited in order to ob- 
tain a trapping potential whose shape depends on the 
internal state of the atoms. As an example, we consider 
an atomic mirror like the one recently realized Q from 
a conventional video tape with sinusoidal magnetization 
M = (Mo cos[kMx], 0, 0) along the a;-axis. The period 
of the pattern, 2n/kM, can be as small as 1/xm with the 
system studied in ||, or even close to 100 nm using ex- 
isting magnetic storage technologies. In order to get a 
microscopic trapping potential pT| , it is necessary to ap- 
ply an external bias field B cxt = (0, By Xt , B° xt ), oriented 
mainly along the z axis, normal to the mirror's surface, 
and with a small component along y in order to prevent 
trap losses due to spin flips occurring at magnetic field 
zeros. In this case the magnetic trapping potential is 

F mF (x) = g Ft i B m F {B^e- 2kMZ cos 2 (k M x) + {B^f 



[B e- kMZ sm(k M x) + B e z xt ] 



y, (37) 



where B = n M (l — er kM& )/2 and S is the tape thick- 
ness. The minima of V mp form a periodic pattern above 
the tape surface, at a height z — ln(^ Af / B ) / k-^ typ- 
ically of the order of some fractions of /im. The spacing 



between two nearest minima along x is just the period 
of the magnetization, 2tt jku- With present-day tech- 
nology, trapping frequencies can range from a few tens 
of kHz up to some MHz. Microscopic electrodes can be 
nanofabricated on the mirror's surface Jl0| , thus allow- 
ing for the design of a potential with the characteristics 
described in Sect. III. 



For the states \a) and \b) we choose the hyperfine struc- 
ture states \a) = \F = 1, mj? = —1) and \b) = \F = 
2,mp = 2) of the 5S1/2 level of 87 Rb, having scattering 



lengths a 



1,1, 



-,ab 



5.1 nm. Several schemes of loading 



atoms into the trap have been envisaged (see for example 
§,|lj). Most of them rely on an intermediate step, where 
atoms can be trapped and cooled without coming in con- 
tact with the magnetic mirror. This pre-loading stage 
can be either a magnetic trap initially displaced from the 
surface, or a different kind of trap (for instance an evanes- 
cent wave mirror, where different internal states can be 
trapped by gravity close to the surface |l3| before the 
atoms are put in the correct states for magnetic trap- 
ping), to be replaced by the electromagnetic microtrap 
with a gradual switch-on of the electric and bias mag- 
netic fields in the final stage of loading |l(| . This could 
also allow for implementing a controlled filling of the trap 
sites by adiabatically turning on the periodic potential, 
in a similar way to that discussed in pTf] . 



B. Results 

1. Time evolution during gate operation 

If both particles are in state |a), there is no interaction- 
induced phase shift, as expressed in Eq. (|27j). The results 
for both particles in state \b) are shown in Fig. 0a, while 
those for differing internal states appear in FigTSb. The 
harmonic potential ensures that the system comes peri- 
odically back to its initial state. In the absence of inter- 
action, the frequency of recurrencies is twice as high for 
IV'fcbW) as i t is for \ipab{t)), as already discussed at the 
III B 2\ . The interaction also makes the two 



end of Sect. 

cases substantially different from each other. Its effect on 
the atomic motion is not dramatic if both particles are in 
state \b): actually, the oscillation period in the presence 
of interaction is increased by only St « 1.4 x 10 _3 T OSC 
with the parameters used here. The collisional phase <pbb 
increases in steps at the times tk = (2k + l)T osc /4, when 
the atoms meet at the center of the well, and remains con- 
stant at intermediate times while they are apart. Note 
that since the particles are indistinguishable the ampli- 
tude for the particles to bounce back during the colllision 
does not harm the perfomance of our scheme. The con- 
tributions of the reflected and the non-reflected part to 
the wavefunction are indistinguishable. What matters is 
whether of not the two-particle spatial distribution ap- 
proaches the initial one, and this is satisfied to a high 
accuracy in our case. 
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The behavior is quite different if the atoms are in dif- 
ferent internal states. The phase shift increases in larger 
steps since the collision is close to the turning point of 
the particle in state |6), near x — xq. Here the velocity 
of the particle is much smaller than at the center of the 
trap and thus the interaction time is longer, allowing a 
larger phase to accumulate. The collision also excites vi- 
brations of the particle in state |a). The resulting loss of 
energy from the particle in state \b) leads to a decreasing 
oscillation amplitude of that particle, and the initial state 
is no longer recovered. This problem can be avoided if 
the potential minimum for state \a) is displaced along the 
transverse direction from the one for state \b) by means 
of an additional electrostatic field [jll) , so that the atoms 
interact if and only if they are both in state \b). 



\O{ip a0 ,t)\ 
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FIG. 2. Dynamics during gate operation: projection of 
the initial state on the state evolved without (left) and with 
interaction (center); interaction-induced phase shift (right). 
Results are shown for different combinations of internal states: 
a) a = P = b; b) a ^ 0. We choose uj = 2tt17.23 kHz 
and oj± = 27rl50 kHz, corresponding to ground-state widths 
a x ~ 82 nm, a± ~ 28 nm, with the initial wells having fre- 
quency u>o = 2uj and displaced by xq = 5a x . Time is in units 
of the oscillation period T osc . 



2. Gate fidelity at T = 

Ideally, the scheme described above should realize the 
mapping 



a a 



-i2tf> a 



WW, 



\a)\b) -» e -*W- + *»+*"»J|o)|6), 
\b)\a) -» e- < ^ + *-+*"»)|6)|a), 
\h\\h\ c-'O^+^iOlMIM 



|15[ to characterize the quality of the gate. F is defined 

as 



F = mm (tr ext { (x\US [\ X ) (xl ® Po] $W \ X ) } 



(39) 



where |x) is an arbitrary internal state of both atoms, 
and \x) is the state resulting from \x) using the map- 
ping ( p8| ) . The trace is taken over properly symmetrized 
motional states, U is the evolution operator for the inter- 
nal states coupled to the external motion (including the 
collision) , S represents symmetrization under particle in- 
terchange and po is the density operator for the initial 
two-particle motional ground state. A straightforward 
calculation gives 



F-- 



1 1-A 2 -B 2 [{1+A 2 )B 2 -4ABC+2C 2 ] cos 2 {<f> bb ) 

2 (1+A){2+B[(1- J 4)B+2C] cos(0 6 i,)}-S 2 (B-C) 2 cos 2 



ibb) 

(40) 



where A = \0(ipH\ r 



B 



|O0M,t)| 2 , c 



lOotybb, T )\ 2 ■ With the parameters quoted above, we 
obtain F ~ 0.99 either by choosing a gate operating time 
t = 7(T osc + 5t), maximizing B, or r = 7T osc , maximiz- 
ing instead A. We prefer this latter choice since, after 
a time r = NT osc = 2Ntt/uj, the j th component of the 
x-wavefunction of an atom in state \b) in the basis of 
eigenstates of Vb(0 < t < r) gets a phase 2N(j + l/2)n 
(here N = 7). This brings some simplifications: e.g., the 
kinematic phases can be written as 

j> a = Nn^— -, cb b = Nn i. 41 

UJ UJ 

The general form of (j>b is much more complicated. Fig. ||a 
shows that after 7 complete oscillations Eq. ([If]) yields a 
phase shift 4>bb{7T osc ) rs 7t due to the interaction, whereas 
the perturbative formula Eq. ( |3l"|) gives 74>b b (T osc ) w 
0.977T. The figure also shows that the overlap \Oo(tpbb, t) \ 
remains close to 1, satisfying Eq. (|l3|). Th e curve has 

71 sig- 



[VB1 



local minima at the times tk defined in Sect, 
nailing that a collision is taking place, and shows a global 
decrease due to the accumulating delay of the interact- 
ing motion with respect to the noninteracting one. The 
fidelity turns out to be 



|Oo(V%6,t)|cos[0»(t)]}. (42) 



(38) 



3. Gate fidelity at T 



where <p a and are the phases due to the time evolution 
without taking into account the interaction. We assume, 
as above, that the trapping potential is designed to pre- 
vent the atoms interacting if they are in different internal 
states. Therefore we set <p a b = in Eq. ( |38| ) and consider 
only 4>bb in the following. We use the minimum fidelity F 



In order to compute the temperature dependence F(T) 
of the fidelity, the density matrix for the motional degrees 
of freedom in Eq. (p9|) has to be replaced by 



Pext (T) = Y,P ln (T)\l) R (l\® 



(43) 



/ . n 
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which coincides with p$ at T = 0. Here we have intro- 
duced the eigenstates \1)r for the center of mass and \n) r 
for the relative motion. The probabilities Pi n (T) for oc- 
cupation of the CM and relative motion excited states 
are calculated assuming for each atom a thermal distri- 
bution corresponding to temperature T, as expressed by 
Eq. @. We obtain 

F(T) = i{l - P m(T) | Oo(^(„),r) | cos [</>(„) (t)]}, 

(44) 



as already discussed in Sect. [II B 3. Consequently the 



l.n 



where 



?=-i,+i 

x H, 



Vn!2«+! 



~2?T 



(2x + ?r) 



(45) 



In particular, ^Vq) = ipbb an d 0(o) = ^fcb- The correspond- 
In) (i) are shown in 



ing interaction-induced phase shifts 
Fig. |b,c. 

<t>{n) (*) 



\OoW in) ,t)\ 




FIG. 3. Dynamics for both atoms in state \b), with rel- 
ative-motion excitations: a) n = 0; b) n = 1; c) n = 2. 
On the left: interaction-induced phase shift; the crosses re- 
fer to the perturbative result from Eq. (po|), explicitly given 
by Eq. for n = 0, and evaluated numerically for n > 0. 
On the right: projection of the evolved state on the corre- 
sponding state evolved without interaction. Trap p aram eters 
have the same values as in Fig. ^| and satisfy Eq. ( B5 ) since 
clqujo I '{Axqw) = 0.07 in this case. 

The discrepancy between the interacting and the non- 
interacting motion increases with n, but nevertheless the 
phase shift </>(„) remains still close to n (Fig. |^b and c), 



fidelity is not rapidly suppressed with temperature. 

For example, one might well be interested in the values 
of F(T) for temperatures up to ksT Hloq. Let us 
therefore define 7 = exp(— Tiwo/fesT) and neglect terms 
of o(j 7 ) in the evaluation of Eq. (jij) to obtain 

1 6 

F(T) « F(0) - - {|OoW<(„),t)| cos |> {n) (T)] 

n=l 

• (ip( n ~i),T)\ cos [^(„-i)(r)]|. (46) 



O, 



This still gives a high fidelity F(T) w 0.96 even at 
fcgT = 2hujQ, for which 7 7 w 0.03. We note that, in 
order to reach such a high fidelity, the timing has to be 
quite precise, with a resolution better than 10 _3 T OSC cor- 
responding to tens of ns in this case. 



V. CONCLUSIONS 

We have shown that entanglement among ultracold 
neutral atoms can be controlled by means of microscopic 
switching potentials. The fidelity for a fundamental two- 
qubit quantum gate turns out to be quite robust with 
respect to temperature: in fact, with the parameters 
quoted below Fig. |, we find F(T) » 0.96 for T rj 3^K 
in the x-motion, while assuming ground-state cooling 
in the transverse directions. We find a gate operation 
time of r « 0.4ms, over which coherence can probably be 
preserved with presently available experimental systems. 
Static microtraps based on available atomic mirrors [p|JTl|l 
provide a good opportunity for a first implementation 
of our scheme. Here nanofabrication technologies allow 
steep potentials to be achieved with small charges and/or 
currents. Trapping fields can be controlled electronically 
in a fast and accurate way JlC[ ] . 

Some problems remain to be addressed. To perform 
even a single gate operation, the trap should be loaded 
with exactly one atom per well. Read-out should be done 
possibly without removing atoms from the trap. In order 
to build up more complex operations, gates should be ar- 
ranged in a periodic structure where coherent atom trans- 
port may take place between different locations. This 
would permit gate operations either on one pair of atoms 
at a time, or on several pairs in parallel, a fact which 
could be exploited for efficient implementation of quan- 
tum error correcting schemes and fault-tolerant quantum 
computing fl6||. This will be the subject of future work. 
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with cu±(t) = y/u) 2 + uj 2 ± (uj 2 — uj 2 ) cos (cut). 

This result for the relative motion should be compared 
to the actual evolution in the presence of interaction, 
which cannot be computed analytically. If the particles 
are in different internal states we also have to resort to 
numerical methods. 



APPENDIX A: TIME EVOLUTION 

1. Analytical calculation 

If both particles are in state \b) we start from the 
Hamiltonian Eq. (p8|), neglect the interaction term, and 
solve the Schrodinger equation. We find (omitting the 
internal state indices bb) 



2. Numerical calculation 

a. Particles in the same internal state 

We write the state vector as a sum over the eigenstates 
|n) of a harmonic oscillator of mass /x and frequency uj, 



(A8) 



lpCM(R,t 

where 

n(t) = ■ 

<f>CM.(R,t) 



irh 



a l<f>CM(R,t) 



(Al) 



UJ 2 UJq 



[uj 2 cos 2 (uit) + ljq sin 2 (ut)] ' 
2 



(A2) 



MQ(i) ujI - uj 2 , . , , , ujt 
— ^-r 1 — R 2 cos sin (cut) - — 



1 

-— arctan 
2 



(luo — uj) cos(ujt) sin(wt) 
ui cos 2 (cot) + u>o sin 2 (cut) 



From Eqs. (26a) and (Al) it follows 

>\2 



|0(VcM,t)| = 



1 



4cu 2 uj 2 



■ sin 2 (u>t) 



(A3) 



(A4) 



If the particles did not interact, the relative motion would 
be 



where 
(f> le i(r,t) ■■ 



Mi" 



i^ Tc i(-r,t)-*^-[r+2x cos{u,t)Y 



Ujt 

2 



(cjq—uj) cos(ujt) sin(ujt) 
uj cos 2 (ijt)-\-LjQ sin 2 (uJT) 



2un(t) . , , 



(A6) 

(^0 — ^ 2 2 2~\ / \ 2 
j r +u) x )cos(uit)+uJ x r 



The overlap between the states Eqs. ( |26b| ) and (A5) is 



v 2„2 ain 2, 
P ■ 



2 cos 


2 / 2 1 2\ 


t 2 




h w 2 (t)w 2 (t) 


Xq 




4mu, 
h 

e 


co3 2 (ct) , sin 2 (^t) 

[ - 2 (t) 1 «iw _ 


w 2 


a: 2 



(A7) 



and approximate the potential by a truncated sum 



k,l 
iV max 

= £ ^(o)^(o)ife)ai, 



(A9) 



where ip n ( x ) — (x\n). We have checked that the final 
result is independent of iV max , with N max of the order of 
some tens. The Schrodinger equation for |i/> ro i(i)) gives 



c n (t) = -i2a?u ± i>* n (0) J2 MOy (n - l)uJt ci(t), (A10) 



1=0 



which we solve numerically for c n (t) w ith < n < A n 
The initial conditions, from Eq. (26b), read 



Cn(0) 



__^o^, z .2 1 . 

e R (^o+^) (wqo;)' (ujq—uj 



Vn!2™ V w o + w V^o + w 



(All) 



2m oJcj 2 a: 2 
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b. Particles in different internal states 

In order to solve the Schrodinger equation for the 
Hamiltonian Eq. (|34|) we decompose the state vector 

llM*)) = E ^ lU+k+1)Qt cMt)\j) R \~k) r , (A12) 

(where now tpj( x ) — ( x \j) are the eigenfunctions of a 
harmonic oscillator with frequency u) and mass m) and 
obtain for the coefficients 







c jk (t)=i 



{c i+1 , fe+1 (t)e- i2a V0' + l) + 



2 (w 2 + w 2 ) 

+ Cj - hk+1 {t)^j {k + 1) + c J -_i, ifc _ 1 (t)e i2fi>t > /i* 

°e 2 c Jfc (t) 



+ c J+ i, fe _i(t)v / (j + l)A: + 



2m£j 



e* 5t (c,_i lfc (t)Vi-q J -,fc-i(t)V^) 
+ e- iS *( Cj - +1)fe (t)V7Tl - ^ >Jk+ i(t)Vfe + l)] } 

-iV2<v^(-e)E^(-^)' 



(A13) 

where £ = x uj 2 / y/2uj 2 . 

This can again be solved numerically, starting from the 
initial conditions, derived from Eqs. (f26a|), (26b), 



Cjfe(O) = 



■ LV 



(A14) 



APPENDIX B: INTERACTION PHASE SHIFT 
FOR EXCITED INITIAL STATES 

Let us consider two bosonic atoms in the same internal 
state | a), but in two different single-particle motional 
states \<p~) and \ip+) with vanishing overlap. The initial 
motional state has the form 



l¥>(0)> = 



+ \^+)\<p-) 

V2 



(Bl) 



We assume that: (i) the particles move against each 
other, come in contact during a certain time interval 
[U,tf] and then separate again; (ii) the velocity of each 
particle and the shape of its wavefunction do not vary 
during the interaction. Thus for U < t <tf we write: 



(xj. | <p_(*)) = ip'(xi - vt), 
(x 2 \ip + (t)) = Lp"{x 2 +vt), 

where v is a positive constant. It follows 

1 f tf 

4> aa (T osc ) w - / dt (ip(t)\u aa (xi : x 2 )\(p(t)) 



(B2a) 
(B2b) 



dt 



+ OC 



dxi 



x -vt)f \cp"( Xl + vt)\ 2 

^± / dxdyW{x)\ 2 W'{y)\ 2 
2a b s b uj 1 _ 



where a change of variables x = x\ — vt, y — xi + vt has 
been introduced, and the limits of integration in t have 
been exten ded t o ±oo since the single-particle wavefunc- 
tions Eqs. ( B2a )-( B2b ) overlap just for a finite time. The 
result turns out to be independent of the initial state. We 
can compare it to Eq. ( pl[) , which was obtained in the 



harmonic potential Eq. (21b) starting from the single- 
particle states \ip±) instead of \<p±). In this case 



dt{il)±\e% nht x e 



X UJ, 



(B4) 



and the atoms collide twice during one oscillation period. 
Therefore the collisional phase Eq. ([5l]) should be twice as 
big as Eq. (B3). This is true provided that the maximum 
velocity for the atomic motion in the well Vb(x, < t < r) 
is large with respect to the analogous quantity for the 
ground-state motion in the wells v a (x,t), i.e. if 



Xqlo > a w /4. 



(B5) 
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